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ABSTRACT 

In earlier work we showed that a dark matter halo with a virial mass of 10^ Mq can survive feedback 
from its own massive stars and form stars for > 100 Myr. We also found that our modelled systems 
were consistent with observations of ultrafaint dwarfs (UFDs), the least massive known galaxies. Very 
metal-poor damped Lyman-a systems (DLAs) recently identified at z ~ 2 may represent the gas 
that formed at least some of the observed stars in UFDs. We compare projected sightlines from our 
simulations to the observed metal-poor DLAs and find that our models can reach the densities of the 
observed sightlines; however the metallicities are inconsistent with the single supernova simulations, 
suggesting enrichment by multiple supernovae. We model two scenarios for the history of these 
systems. The first explains the gas abundances in DLAs by a single burst of star formation. This 
model can produce the observed DLA abundances, but does not provide an explanation as to why the 
DLAs show suppressed [a/Fe] compared to the stellar population of UFDs. The second scenario splits 
the DLAs into a population which is enriched by a single burst, and a population that is enriched by 
a second burst after the accretion of metal-poor gas. In this scenario, the suppressed average [a/Fe] 
in DLAs compared to UFDs results from enrichment of second-burst systems by Type la supernovae. 
Subject headings: first stars, reionization, galaxies: abundances, galaxies: dwarf, galaxies: formation, 
stars: Population H 


1. INTRODUCTION 
The very metal-poor Damped Lyman-a s ystems 


(DLAs) with [Fe/H] < —2 recently discovered by Cooke 
et al. (2011a|b 2013 20141 have chemical signatures sug- 


gestmg that they have experienced only a few enrichment 
events, making them promising candidates for probing 
the chemical signatures of the first generations of stars. 
DLAs absorb light from a bright background source, al¬ 
lowing the detection of gas that is too faint to observe in 
emission. The most distant and metal-poor DLAs may 
be the environments of the formation of some of the first 


logical hydrodynamical simulations (iPontzen et al.|2008 

Fumagalli et al. 

2011 

Cen 

2012 

Bird et al. 

2013[), the 


(Krogager et al.p012 Jorgenson fc Wolf^|2014) suggest 


that the host galaxies of higher metaliicity [Z 
D LAs are most l y 10^ — 10^^ Mq dwarf galaxies. 


1/50) 


Cooke et al. (2014) found that the kinematics of the 
very metal-poor DLAs agreed with those of Milky Way 
dwarf spheroidal galaxies. Furthermore, the decline in 
[a/Fe] for [Fe/H] > -2.0, which is usually a ssumed to indi - 


cate enrichment from Type la supernovae (Tinsley 1979), 
is similar to that observed in the dwarf spheroidals, in- 
dicative of a similar star formation history. This differs 
from the Milky Way halo, which shows approximately 
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constant [a/ Fe] for —3.75 < [Fe/H] < —0.75. Cooke 


et al. (2014) therefore suggest that the very metal-poor 


DLAs correspond to the least massive systems that can 
form stars at z ~ 3. 

Over the past decade, a number of galaxies with lumi¬ 
nosities below 10® Lq have been discovered in the local 
universe. Known as ultrafaint dwarfs (UFDs), these sys¬ 
tems are promising candidates for probing the chemical 
signa tures of the first generat ions o f stars 


I'rebel et al.||2014|). The majority of Ul'Ds contain only 

ancient stars, with the spread of ages being < 2 Gyr 

(IMuhoz et al. 

2010 Brown et al. 2012 Vargas et al. 12013 ). 

Brown et al. 

(201H found that the stellar populations 

ot five Ul'Ds 
1 and Coma 
that a global 
star formatio 
galaxies with 

: Hercules, Leo IV, Ursa Major 1, Bootes 
Berenices showed similar ages, suggesting 
event such as reionization truncated their 

n. Weisz et al. (2014b) found that most 

M* < lO'^ Mq formed > 80% of their stars 


before z ~ 2, although there was significant variation be¬ 
tween individual galaxies. For example. Canes Venatici 
H formed stars until ~ 8 — 11 Gyr ago, while Hercules 
and Leo IV formed 90 % of their stellar mass by 11- 
12 Gyr. The star formation history of UFDs indicates 
that at least some of them were forming stars at z ^ 3 
and therefore must have contained neutral gas at this 
time. This suggests the possibility that very metal-poor 
DLAs trace gas in UFDs at the time they were forming 
stars. 
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With current instruments, it is difficult to observe 
UFDs beyond the local Milky Way + M31 system, al- 
though there has b een one detection in the Virgo cluster 
(Jang & Lee||2014). The relationship between low mass 
tossil galaxies that formed in the early universe and the 
observed UFDs today is therefore complicated by the in¬ 
fluence of environment. Current observations cannot dis¬ 
tinguish between systems which are true fossils of reion¬ 
ization and those that had their star formati on quenched 


by oth er processes, such as tidal stripping (Weisz et al. 
2014a). The UFDs that may be traced by DLAs need not 


be close to a host galaxy and could provide an opportu¬ 
nity to study the smallest galaxies in a less complicated 
e nvironment. 

Weisz et al. ( j201^ studied 13 galaxies identified by 


Bovill & Ricotti 1 2011) as fossil candidates, finding that 


Hercules and Leo iV were the only strong candidates 
based on their star formation history. The Cooke DLAs 
trace gas at z ^ 2 — 4, indicating that the systems they 
trace were not permanent ly quenched by reio nization. 
This is consistent with the Weisz et al. (2014b) star for¬ 
mation histories for systerns such as Canes Venatici II, 
for which the best fit model formed the bulk of its stars 
at z ~ 2. The mass of Cane s Venatici II is 1 .4 x 10® Mq 
within the half-light radius ( Wolf et al.|2010 ), suggesting 
that even very low-mass systems can form stars at the 
redshifts of the Cooke DLAs. 

In this work we use the hvdr qdvnamical simulations 


of a M yi,. = 10^ Mq galaxy from Bland-Hawthorn et al. 
(2015 ), in conjunction wit h the chemical evolution model 


from 


Webster et al. (2014), to investigate possible scenar¬ 


ios for the history and evolution of DLAs. Two main pro¬ 
cesses prevent galaxies with lower masses from forming 
or surviving in the early universe. The first is the en- 


Hawthorn et al. 

2011 

2015 

Webster et al.[ 


formation masses of 10^ Mq and lower, showing that 
the limit for stars to retain gas and form stars in the 
face of feedback from their own star formation is Mvir ~ 
10®'® — 10^ Mq. Furthermore, our models were able to 
reproduce most features of the observed relationship be¬ 
tween [a/Fe] and [Fe/H] in ultrafaint dwarf galaxies. 

The second process that can remove neutral gas from 
galaxies is the epoch of reionization, when the radiation 


in the universe 

. It is often claimed (Reesll 19861 Barkana 

& Loeb 

1999 

Gnedin 

2000 

Ukamoro et al.| 2008p that 


neutral gas is unable to cool onto halos -wit h masses 


< 10®Mq. H owever, more recent simulations ( [Bovill fc 
Ricotti|[2009 ) suggest that this is not a hard limit, with 


the influence of environment meaning that at least some 
lower mass systems can survive the epoch of reionization. 
Ricotti & Gnedin|(|2005|) used cosmological simulations to 
predict the existence or an undetected population of ul¬ 
trafaint galaxies, around the same time as the first UFDs 
were observed. 

We now consider whether 10^ Mq halos can explain 
the number of observed very metal-poor DLAs. For each 
10^^ Mq halo, there are ~ 10^ halos with M < 10^ Mq. 
At z ~ 2 — 4, the virial radi us of a 10^^ Mq halo wa s 
^ 120 kpc (calculated from Bryan & Norman 1998), 
while the virial radius of a 10' Mq halo was ~ 0.5 kpc 


(see Figure 1 and Bland-Hawthorn et al. (2015|)). The 
ratio of cross-sectional area is proportional to tne ratio 
of their radii squared, which is ^ 2 x 10“®. However, 
10^ Mq halos are 10"* times more common, so the total 
cross-sectional area is ^ 20% that of the 10^^ Mq halo. 
While this is only an approximate order of magnitude ar¬ 
gument, it suggests that DLAs from 10^ Mq halos should 
be frequent enough to be observed. 

In Section we will summarise the simulations, al- 
thor^h we refer to our previous work for the details. Sec- 
tionj^compares the Cooke DLAs to our models. Possible 
scenarios for the history of the DLAs and their link to 
present-day UFDs are discussed in Section followed by 
our conclusions in Section 

2. SIMULATIONS 

The simulations used in this work a re described in de- 
tail in iBland-Hawthorn et al. (2015|) fc Webster et al. 


(2014).^ In [Bland-Hawthorn et al. ( 2U15[) , we pre¬ 


sented high-resolution simulations performed using the 
hydro/ionization code Fyris Alpha showing the effect of 
a 25 Mq star on low-mass dark matter halos. An ionised 
region is created around the star prior to the supernova, 
resulting in the supernova having a much greater effect 
than for a lower mass star. The thermal and ionisation 
structures were calculated with the MAPPINGS IV ioni- 


sation code, with the AT LAS9 atmospheric g r id ( ^astelli 
& Kurucz 2004) and the Meynet & Maeder (2005) evo- 


lutionary tracks used to model the proge mtbr star. The 
full de tails can be found in Section 2 of lWebster et al.l 

([Mij. 

As a check for convergence, we investigated the amount 
of energy lost given a single-level resolution model, a two- 
level model with three times more cells in the central 
region, and a three-level model with the resolution again 
higher by a factor of 3. If the resolution is insufficient, the 
gas has a tendency to overcool. The three models showed 
only minor differences, and we therefore used the two- 
level model with the higher resolution level of 216® cells 
having a resolution of ^ 2 pc/cell. The other possible 
effect of insufficient resolution is that the metals will mix 
too efficiently. 

The inclusion of radiative cooling and an inhomoge¬ 
neous interstellar medium reduces the coupling between 
the supernova energy and the gas, such that halos with 
masses Mvir = 10®'®“^ Mq retain a large proportion of 
their baryons in the face of the ionization and supernova 
of the 25 Mq star. The location of the star in the halo is 
also important, with supernovae occurring away from the 
centre resulting in lower levels of enrichment and greater 
retention of dense gas. 

The simulations start with an intergalactic medium 
enriched to [Fe/H] = —4 by the first stars. Such low 
[Fe/H] has been obser ved in the IGM at z ^ 3 — 3.5 
(Fumagalli et al. [2011 ). Our enhanced C, O and a el- 
ements ([a/Fe]init = 0.7) give Zjnit = 10“®'^ Zq. We 
take our initial metallicity to be the critical threshold for 
low-mass star formation, such that low mass star forma¬ 
tion proceeds only after the initial 25 Mq star enriches 
the gas. While our thresho ld fits the critical me t allicity 
Zciit 10“^ — 10“® Zq of Smith & Sigurdsson[ (2007), 
other authors have suggested t hat this threshold couM 
be as low at Zcrit ~ 10“® ^ (Schneider et al. 2006). 


As discussed in Webster et al. (2014), removing the as- 
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sumption that low-mass stars do not form at our starting 
metallicity does not have a large impact on the results, 
as it affects only the first 10 Myr of the 600 Myr star 
formation history. 

The gas mass within the scale radius of the Einasto 
dark matter potential is set to be 10% of the dark 
matter mass within this radius. The overall baryon 
fractions Mgas.vir/Tfvir are « 12%. The model of an 
Mvir = 10^ M© halo contains 2.34x10® Mq of gas within 

Ts and 1.5x10® Mq within Cvir- _ 

In the hydrodynamic al simulations presented in|Blancb| 


Hawthorn et al. (20151, the evolution of the gas after a 
single supernova was traced for 60 Myr. Longer time pe¬ 
riods were not possible in these high-resolution simula¬ 
tions due to the effect of boundary conditions, while dif¬ 
ficulties involved with studying the interaction between 
su pernovae restricted th e models to a single supernova. 
In Webster et al. (2014), the density and metallicity of 
the gas after a single supernova was used as a template 
for the condition of the gas after subsequent supernovae. 
This allowed the density and metallicity of the gas to be 
tracked for 600 Myr. We model our halos in isolation, 
so do not include infall of diluting gas from the IGM, 
although we do consider the effects of dilution in Sec¬ 
tion 4.2, where the second burst of stars forms after the 
accretion of metal-poor gas. 

The hydrodynamical simulation does not trace molecu¬ 
lar cooling as would be required for the gas to cool to star 
formation temperatures. Instead, stars are allowed t o 
form in the gas using the method of Argast et al. (2000), 
in which 10^ cells are randomly selected, with each given 
a probability of forming a star proportional to the square 
of the density of the gas in the cell. Th e masses of the 
stars formed were selected by sampling a Kroupa (2001) 
initial mass function. After each supernova, the density 
distribution of the gas is reset to the distribution just af¬ 
ter the supernova in the single-supernova hydrodynam¬ 
ical simulation and the density then evolves as for the 
hydrodynamical simulation. The metallicity is treated 
in the same way, except that it is added to rather than 
reset, such that the enrichment from each supernova adds 
to the enrichment from all previous supernovae. 

The yields used to deter mine [Fe/H] & [a/Fe] are from 


Woosley & Weaver (1995), interpolated and extended to 
8 Mq , taken to be tne lowest mass star that ends its life 
as a supernova. After 100 Myr, Type la supernovae occur 
with a probability e quivalent to th e Type la supernovae 
rate determined by Jimenez et al. (2014), scaled to the 
star for mation rates of our sys tems. 'Lhe Type la yields 
are from Iwamoto et al. (1999). 

Fig. shows the halos used in this work. The Mvir = 
10®'® Mq, 10® Mq, 10®'® Mq and 10^ Mq models are 
referred to as M55, M60, M65 and M70 respectively. The 
scale radius, which contains nearly all the gas sufficiently 
dense to be observed as DLAs, ranges from 33 pc for M55 
to 151 pc for M70, while the virial radii range from 200 pc 
for M55 to 630 pc for M70. These halo properties are at 
z = 10. As shown in Figure]^ halos with Mvir = lO®”"^ at 
z = 10 grow by an order of magnitude in both mass and 
radius by z = 3. Their virial radii grow from 0.3-0.6 kpc 
to 1-3 kpc. The increase in halo mass will improve the gas 
retention compared to that of the first supernova. Our 
assumption that the first supernova occurs at z = 10 



X (pc) 


Fig. 1.— The physical sizes of the modelled halos at z = 10. M55 
denotes a dark matter virial mass of 10^'^ Mq, MOO has Mvir = 
10^ Mq etc. The inner circles represent the scale radii, while the 
outer circles are the virial radii. For the off-centred explosion, the 
supernova is placed at the radius enclosing half the gas inside the 
scale radius, which is a pproximately halfway between the centre 
and the scale radius. See |Bland-Hawthorn et al.| (|2015|l for the full 
set of parameters. 


Redshift (z) 

20 10 5 2 1 0 



Fig. 2.— Evolution of the virial mass and radius for a halo with 
Mvir = 10^ Mq halo at 2 = 10 (indicated by the dotted line) from 
2 = 20 to the present time. This plot was produced using dat a 
from 5000 runs of the tree merger code of|Parkinson et al.| (|2008||. 
The dashed lines show the range which 6y7o ot halos tali within. 


is supported by Power et al. (2014), which shows that 
halos with present-day masses Mvir < 10® Mq reach the 
atomic cooling threshold at approximately this redshift. 

The suffixes we use in the rest of this work represent the 
type of model: CCH represents a model with a central 
explosion, a clumpy ISM and a preionisation phase, OGH 
is the same except with an off-centred explosion, while 
GGC and OCC represent models without a preionisation 
phase. For example, the M650GG model is a halo with 
Mvir = 10®'® Mq, which experiences an off-centred explo¬ 
sion (defined as a supernova at the radius enclosing half 
the gas mass within the scale radius) in a clumpy ISM, 
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Fig. 3.— The distribution of H/ column density vs projected 
metallicity 25 Myr after the explosion of a 25 Mq star for the 
models with a preionization phase. The vertical dashed line in- 
dicates the definiti on of the minimum DLA column density. The 
|Cooke et al.| )|2014[l DLAs are also plotted with their error boxes. 


Fig. 4.— The distribution of H/ column density vs projected 
metallicity 25 Myr after the explosion of a 25 Mq star for the 
models without a preionization phase. The vertical dashed line 
indicates the defini tion of the minimum DLA column density. The 
|Cooke et al.||^014[l DLAs are also plotted with their error boxes. 


with no prei onisation phase. Only the Bland-Hawthorn) 
et al. (20151 models with a clumpy ISM and radiative 
cooling are used, as these are the most realistic. 

3. MODELLED DLAS 

In this section, we discuss the projected DLAs from the 
simulations discussed briefly in Section and in detail in 
Bland-Hawthorn et al. (2015). The column densities and 
metallicities result troin enrichment by a single supernova 
from [Ee/H] = —4. The metal-rich supernova ejecta were 
traced with a scalar variable, advected passively, such 
that the local enrichment of the gas could be determined 
at each time-step. Figs. 3 5 show the distribution of col¬ 
umn densities and met^icities along projected lines of 
sight, while Figs. |6][^show spatial maps of column den¬ 
sity and metallicity. The state of the gas after 25 Myr is 
used, as this is sufficient time for th e gas to recover from 
the supernova ( Webster et al.| 2014). 

In this section, we assume that our halos exist in isola¬ 
tion and do not consider the complex possibility of merg¬ 
ers between low-mass gas-rich halos. A merger could in¬ 
troduce a fresh supply of low-metallicity gas, resulting 
in dilution as in the scenario considered in Section 4.2. 
This would result in higher column densities and lower 
metallicities, such that more sightlines would fit the DLA 
threshold. This would make it less likely that the sys¬ 
tems could reach the metallicities of the DLAs with only 
a single supernova. The merger could also induce a burst 
of star formation. The effect of this is less certain, with 
possibilities including little or no change to the results, 
a more rapid increase in metallicity, or the removal of a 
large proportion of the gas in the system due to multiple 
supernovae close together in time. 

We first discuss the models with a strong preionisation 
phase from a 25 Mq star, then the models with a weak 
preionisation phase, consistent with an M < 15 Mq star. 
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Fig. 5.— The distribution of H/ column density vs projected 
metallicity 25 Myr after the explosion of a 25 M 0 star for M60 with 
a preionisation phase and M55 without a preionisation phase. The 
vertical dashed line i ndicates the definitio n of the minimum DLA 
column density. The |Cooke et al.| ||2014|| DLAs are also plotted 
with their error boxes. 

3 . 1 . Strong preionisation phase 

As can be seen in the top panel of Fig.[^ the M70CCH 
model retains a large amount of dense gas within the 
scale radius (150 pc), with the entire region having a 
column density Nhj > 10^®'^®cm“^. A number of sight¬ 
lines within 50 pc of the centre are dense enough to be 
classified as DLAs. The metallicity of the densest region 
is [Fe/H] ^ —3.5, with some gas reaching [Fe/H] = —3.0. 
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The bottom-left panel of Fig. shows that the highest 
density regions are the most metal-rich. This results from 
the central location of the supernova, such that the met¬ 
als are initially deposited into the dense central region, 
along with it being the most massive model, such that 
more metals are retained close to the centre rather than 
reaching the halo or escaping into the IGM. No other 
model shows this trend. Several sightlines in M70CCH 
reach column densities greater than 10^^ cm“^ 

The M70OCH model is shown in the bottom panel of 
Fig.[§ Unlike in the central case, the dense regions do 
not coincide with the higher metallicity regions. The 
metals from the off-centred explosion do not reach the 
central region and the densest gas is therefore not en¬ 
riched. The bottom-right panel of Fig.|^show that there 
are a few sightlines outside the scale radius that reach 
[Fe/H] = —2.5 with column densities close to the mini¬ 
mum required to be defined as a DLA. These metallicities 
are higher than is seen in the central explosion model, but 
this is because there is less hydrogen gas in this region. 
The metal escape fraction is actually higher than in the 
central case because the metals can escape through the 
lower density regions away from the centre. Most of the 
dense gas has [Fe/H] < —3.5 and the sightlines with col¬ 
umn densities greater than 10^^ cm“^ show almost no 
enrichment. 

The M65CCH model shown in the top panel of Fig. 
has lower column densities than M70CCH with no sight¬ 
line meeting the definition of a DLA. This model reaches 
higher [Fe/H] because there is less hydrogen. The to¬ 
tal amount of metals is less than in the M70 case, as 
a greater proportion of the metals escape into the IGM 
as a result of the lower halo potential. The metallicity 
distribution in the top-left panel of Fig. [^ does not show 
an increase in metallicity with increasing column density 
because the enriched regions are those most affected by 
the supernova and therefore lose a significant amount of 
their dense gas. 

The M650GH model in the bottom panel of Fig.j^and 
the top-right panel of Fig.^shows similar features to the 
M70OGH model, although the column densities are lower 
and the metallicities are higher for the same reasons as 
discussed for M65GGH. A few sightlines exceed the DLA 
threshold, however none of these are enriched above the 
starting metallicity [Fe/H] = — 4. The lines of sight 

with the highest metallicity [Fe/H] ~ — 2 have column 
densities of ® cm“^. 


3.2. Weak preionisation phase 

The M65GGG model is shown in the top panel of Fig. [8 
The column density distribution is similar to M 70GGH 
model, which is consi stent with our finding in [Bland- j 
Hawthorn et al. (2015|) that a decrease of 0.5 dex in halo 
mass has a similar enect on gas retention as switching 
off the preionisation phase. [Fe/H] is higher than in the 
M70CCH model because the gas is less dense while the 
chemical yields from the star are the same. The metallic¬ 
ity distribution in the bottom-left panel of Fig. [^ shows 
almost no variation with column density. 

The M650GG model in the bottom panel of Fig. [^ 
also shows a similar column density distribution to its 
M70 counterpart with preionisation, but once again the 
metallicity distribution is different. The lower gas den- 
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Fig. 6.— M70 model 25 Myr after the explosion of a 25 Mq 
star. Top: Central explosion. Bottom: Off-centred explosion. The 
projected column densities are shown in the contour maps on the 
right, while the metallicities are shown on the left. 


sities in the M65 halo allow some metals to reach and 
enrich the dense gas in the centre. This can be seen in 
the bottom-right panel of Fig. where some of the gas 
within the scale radius reaches^e/H] > —3. 

The lack of a preionisation phase means that even 
M60 halos can retain sufficient gas to reach the DLA 
threshold, as is shown in Fig. [^ In M60CCG, a few 
sightlines reach column densities of cm“^ with 

[Fe/H] = — 3.0. Unlike in the M65 case, there is no 

gas within the scale radius with [Fe/H] < —3.5, showing 
that the supernova has enriched all the central gas. As 
shown in the top-left panel of Fig. this is the only sys¬ 
tem other than M70GGH for which metallicity increases 
with column density. 

The M60OCC model shows a low density, high metal¬ 
licity region near the supernova which is a feature of all 
the models with an off-centred explosion. The top-right 
panel of Fig. [^ shows that the regions with column den¬ 
sities exceeding the DLA threshold remain unenriched. 

Fig.[^shows the lowest mass systems we modelled. The 
M60CII models and the M55GG models do not retain 
dense gas in the face of the supernova explosion, with 
the maximum column density being cm“^. In 

the off-centred cases the densest gas remains unenriched, 
while the models with a central explosion retain gas with 
a column density of 10^® cm“^ with [Fe/H] ~ — 2.5. 


3.3. Comparison with the Cooke DLAs 

In this section we compare our models above to the 
23 observed very metal-poor DLAs presented by [Gooke 


''e/H 


2014). These DLAs have redshifts of z = 2.1—4.5, 
ranging from —3.48 to —1.86, and column densi- 
= 19.6-21.4. Figs. [|[5 
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Fig. 7. — M65 model 25 Myr after the explosion of a 25 Mq 
star. Top: Central explosion. Bottom: Off-centred explosion. The 
projected column densities are shown in the contour maps on the 
right, while the metallicities are shown on the left. 
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Fig. 8.— M65 model 25 Myr after the explosion of a star with¬ 
out a preionization phase. Top: Central explosion. Bottom: Off- 
centred explosion. The projected column densities are shown in 
the contour maps on the right, while the metallicities are shown 
on the left. 
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Fig. 9.— M65 model 25 Myr after the explosion of a star with¬ 
out a preionization phase. Top: Central explosion. Bottom: Off- 
centred explosion. The projected column densities are shown in 
the contour maps on the right, while the metallicities are shown 
on the left. 


modelled DLAs, with the highest metallicities in our 
models resulting from low amounts of neutral hydrogen 
and therefore low column densities. This suggests that 
enrichment from a single supernova is generally not suf¬ 
ficient to explain the observed DLAs. We will discuss 
extended star formation in Section However, a few 
DLAs could be explained by enrichment from a single 
star. 

Two systems show [Fe/H] ~ —3.5. The gas masses are 
not known for these systems, although one has an upper 
limit of Mwnm < 6.3 x 10® M©. In our M70CH mod¬ 
els, which have a gas mass of 10® Mq within and 
2 X 10® Mq within r^, the 25 Mq star enriches the sur¬ 
rounding gas to metallicities ranging from [Fe/H] = —4.0 
to —3.0. The two lowest metallicity DLAs are therefore 
consistent with enrichment by a single supernova in a 
10^ Mq halo. The M65CCC model also shows sightlines 
with the correct column density and metallicity. These 
two DLAs may therefore trace gas polluted only by a 
single supernova. This level of [Fe/H] is con sistent with 
that suggested by Bromm & Yoshida (20111 for the first 
galaxies, so they may even only be enriched by Popula¬ 
tion HI stars. 

Two other DLAs could be explained by enrichment 
from a single supernova if they are lower mass systems. 
HS0105-I-1619 and J2155-I-1358 have column densities 
« 10^® ® cm“^ and [Fe/H] « —2.1. This is consistent 
with the M65CH models. However, both show [a/Fe] 
<0.3, which is difficult to explain given enrichment from 
a 25 Mq star. The chemical abundances are more con¬ 
sistent with the higher column density systems, so it is 
more likely that they have experienced an extended pe¬ 
riod of star formation followed by losing most of their 
gas. The column densities are consistent with the sys- 
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TABLE 1 

Properties of DLAs compared to our model. The masses, 

RADII AND densities FOR THE MODELS ARE TAKEN AT THE SCALE 
RADIUS OF THE ElNASTO POTENTIAL. 


terns that form no further stars in our models. 

The remainder of the DLAs have metallicities to o high 
to be explained b y a single supernova. However, in |Web-| 
ster et al. (2014) we showed that the gas recovers from 
the impact of a 25 Mq star in less than 25 Myr for an 
M70 halo regardless of supernova location. Furthermore, 
the M65 models without a preionisation phase show little 
evidence of disruption or gas loss, suggesting that stars 
less massive than 25 Mq are unlikely to blow out a sig¬ 
nificant proportion of the gas in an M70 system. Our 
DLA models are therefore likely to remain relevant for 
s ubsequent supernova e. 


Cooke et al. (2014) determined the thermal contribu- 
tion to the line broadening for 12 clouds in 9 of the 23 
systems, allowing gas properties to be calculated. The 
total warm neutral medium gas masses range from less 
than 3 x 10^ to 2 x 10^ M©. A comparison between the 
properties of the gas in our model and the properties 
determined for these DLAs is shown in Table 1. The 
neutral hydrogen density n{H), the cloud radius th^ and 
the warm neutral gas mass Mwnm for our model are at 
the scale radius. It is possible that the gas observed in 
DLAs extends beyond these radii, which would result in 
lower densities and higher gas masses. The temperature 
and sound speed are slightly lower in the models than 
for the observations. Overall, the DLAs for which these 
physical quantities could be determined are consistent 
with halos with 10 ^ Mq or slightly higher. 

4. ENRICHMENT OF THE DLAS 

In this section we will discuss two possible enrichment 
histories for the DLAs that show evidence of extended 
star formation. The first scenario suggests that they have 
been enriched only by a single burst of stars, with the 
alternative being that they are a mix of systems which 
have been enriched by one burst and those which have 
been enriched by two bursts foll owing the accre t ion of 


m etal-poor gas. As discu s sed in Webster et al. (2014) 


& Bland-Hawthorn et al. 

(2015 

), only the M70CH and 

the M65CC models retain 

iense 

enough gas for extended 


low mass star formation. This section focusses on the 
M70CH model, usin g the chemical evolution model of 


Webster et al. (2014). 


4.1. Single burst 

The single burst scenario explains the relationship be¬ 
tween [a/Fe] and [Fe/H] in the following way: 


1. Gas enriched to [Fe/H] 
matter halo. 


-4 condenses onto a dark 


crease in [Fe/H] and a decrease in [a/Fe]. [a/Fe] 
declines at low [Fe/H] due to Type 11 supernovae 
from lower mass (8-1 5 Mr?)) stars, which have 



M70 

M65 

< DLAs > 

DLASmax 

DLASmin 

Tsas (K) 

4860 

2300 

9600 

17000 

5600 

log(n(iT)/cm ®) 

-0.13 

-0.07 

-1.0 

-0.35 

-1.3 

r/Li (pc) 

150 

91 

220 

1270 

32 

ATwnm (IO^Mq) 

2.3 

0.7 

2.5 

220 

< 0.4 

Cs 

5.8 

4.0 

8.0 

10.7 

6.1 


a mean [g/Fe] of 0 .2 (Woosley & Weaver 1995 
Nomoto et al.||2006[ extrapolated to 8 Mq) This 
decline commences at lower [Fe/H] than in our 
model. 


3. After t ^ 100 Myr, Type la supernovae enrich the 
gas, resulting in a decline in [a/Fe] with increas¬ 
ing [Fe/H]. Thi s decline is observed in DLAs with 
[Fe/H] > -2.0 ( jVladilo et aL][MT] ). 


X 



Time (Gyr) 

Fig. 10. — The gas and metallicity evolution as a function of time 
for one realisation of the single-burst scenario in our simulations. 
Top: The total star formation since t = 0. Bottom: The mean 
metallicity of the gas within the (200 pc)® main simulation level as 
a function of time. 

The stellar mass and gas metallicity evolution in one 


2. Low mass star formation commences and Type H 
supernovae enriches the gas, resulting in an in¬ 


simulation run for this scenario is shown in Figure 10 
Our model is stochastic, such that different simulation 
runs will provide different results, but the overall features 
are similar. The star formation rate is on average con¬ 
stant with time, but is bursty, with several short periods 
of little or no star formation following each supernova 
explosion. Each supernova causes a large jump in the 
mean metallicity in the central region, which is caused 
by the deposition of new metals onto the grid, along with 
the lower metallicity gas being pushed outwards beyond 
the inner (200 pc)^ region considered here. The mean 
metallicity in this region then decreases as gas returns 
to the centre, although some of the metals have mixed 
with the gas and the overall metallicity therefore remains 
higher t han before the supernova. 

Figu re 11 shows [a/Fe] vs [Fe/H] for the Cooke et al. 
(2014) DLA observations compared to the gas m our sim- 
ulation, as w e ll as o bservations of stars in UFDs from 
Vargas et al. (2013). The three are reasonably consis- 
tent for [I'e/HJ > —2.5. However, at lower metallici¬ 
ties both our model and the UFD observations show a 
gradual decline in [a/Fe] with [Fe/H], while the DLAs 
show a rapid decline between [Fe/H] ^ —3.5 and —2.8, 
with [a/Fe] then remaining constant until [Fe/H] = —2.0. 
The mean [a/Fe] abundance of 0.25 for the DLAs for 
—3 < [Fe/II]< —2.5 is also lower than the average of 
0.35 for Type H supernovae in a Kroupa or Salpeter IMF, 
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which is observed in sta rs in the halo of the Milky Way 


(Frebel & Bromm 

2012 

)• 

Cooke et al. ( 

201 

4) also notes the suppression in [a/Fe] 


for the ULAs compared to stars in the overall popu¬ 
lation of dSph galaxies. They suggest that this could 
be caused by small number statistics, modelling tech- 
nique s, or physical differenc es between the two popula¬ 
tions. Karlsson et al. (2012) suggested two classes of ex- 


planations tor suppressed [a/Fe] in a star cluster in Sex¬ 
tans. The first explains the low [a/Fe] by contribution 
from Type la supernovae at low [Fe/Hj. This could result 
either from the accretion of low metallicity gas, reducing 
[Fe/H] at a time when Type la supernovae contribute 
to the enrichment, or from a very low star formation 
rate. The seco nd class of explana tions involve only Type 
II supernovae. Aoki et al. (2009) note that a truncated 
IMF at high masses would result in lower [a/ Fe], because 
higher mass st ars yield more alpha elements. [Weidner fc| 


Kroupa (2005) suggest that such an IMF is required for 
systems with very low star formation, which are likely 
to have low cluster masses, limiting the maximum mass 
for a star that can form in the cluster. The final al¬ 
ternative discussed is hypernovae, which produce more 
Fe than normal Type II supern ovae. However, bas ed on 
iron-peak element abundances, Cooke et al. (2013) sug¬ 
gest that very metal-poor DLAs were enriched by stars 
that exploded as core collapse supernovae which released 
1 .2x10^^ erg of energy. 

The question of whether Type la supernovae are re¬ 
sponsible for the suppression of [a/Fe] in DLAs is not 
easily resolved. However, assuming the DLAs are first 
burst systems and a reasonably homogeneous population, 
explanations for the decline in [a/Fe] between [Fe/H] = 
—3.5 and —2.5 involving Type la supernovae are unable 
explain the constant [a/Fe] for —2.5 < [Fe/H] < —2.0. It 
is more likely that for a single burst. Type la supernovae 
start occurring at [Fe/H] « —2.0 as was concluded by 
Cooke et al. ( 2014[ ). 


Given that enrichment by Type la supernovae is un¬ 
likely to explain the low [a/Fe] for —2.5 < [Fe/H] < 2.0, 
the single burst scenario requires that supernovae with 
lower mass yields have enriched the gas observed in 
DLAs. One possibility for this is that DLAs are enriched 
by low-mass (8 — 15 Mq) Type H supernovae to a level 
lower than the average [a/Fe] for a Kroupa or Salpeter 
IMF. One possible explanation for this is that for the rea¬ 
sons discussed above, the IMF is truncated at the high 
mass end. We model an IMF truncated at 20 Mq, with 
the results shown in the top panel of Fig. [^ [a/Fe] re¬ 
mains slightly enhanced compared to the [Fe/H] ~ — 2.5 
DLAs. Reducing the highest mass further would ht the 
lowest metallicity DLAs, but could not simultaneously 
explain the DLAs with [Fe/H] > —2.5. 

Another possible exp l anatio n is a selection effect. In 


Bland-Hawthorn et al. (2015) we studied the effect of 
a 25 Mq star on halos with masses 10® — 10^ Mq. It 
was found that a strong preionization phase evacuates 
the region close to the star, resulting in the subsequent 
supernova having a much larger effect on the ISM than 
it would in the case of a < 15 Mq star with a much 
weaker preionization phase. In particular, a 25 M r^ star 


in a 10®-® Mf7^ h alo can terminate star formation (Web¬ 


ster et al. 2014). It is therefore possible that systems 


which have experienced lower mass stars and therefore 
lower [a/Fe] are more likely to contain dense gas and be 
observed as DLAs. This would have a similar effect to 
the truncated IMF as in the top panel of Figure |T^ 

The above explanations provide a number of possi- 
ble ways that Type H supernovae could result in the 
suppressed [a/Fe] seen in very metal-poor DLAs com¬ 
pared to the Kroupa IMF. However, it does not ex¬ 
plain why the DLAs show lower [a/Fe] compared to 
stars in UFDs, showing [a/Fe] lower by 0.25 dex for 
—3.0 < [Fe/H] < —2.5. None of the observed UFDs 
has mean [a/Fe] < 0.4 for [Fe/H] < —2.5. We therefore 
conclude that if all DLAs are systems forming their first 
burst of stars, either they are not representative of the 
UFD population, or there are systematic differences that 
arise in the process of measuring [a/Fe] between gas in 
DLAs and stars in UFDs. It is possible that as a result 
of evolution, the DLAs are not representative of UFDs as 
observed today. Many of the DLAs contain sufficiently 
dense gas to continue forming stars. In the next section 
we consider a scenario where the difference between the 
two classes of systems is explained by the DLAs being 
observed at different points in their evolution. 


4.2. Two bursts 

We now investigate a scenario in which the DLAs are 
a mix of systems that have formed only one burst of 
stars, and systems that have formed or are forming their 
seco nd burst of s tars af ter the accretion of metal-poor 
gas. Weisz et al. (2014b) investigated the star formation 
history of 38 Local G roup dwarf ga l axies, including 4 of 
the UFDs studied in [Vargas et al.j (|2013). They found 
that Hercules was consistent with having formed all its 
stars in one burst before and during reionization. Leo IV 
also shows only one burst, but likely formed stars until 
z ^ 2. Leo T and Canes Venatici H have a more extended 
star formation history, although Leo T shows no star 
formation from z =1 — 5. The best fit model to Canes 
Venatici H shows a burst of star formation commencing 
at 2 ^ 2.5, while the Icr uncertainty is consistent with 
80% of the stars being formed in a single burst ending 
at z ^ 2.5. This suggests that star formation occurred 
in some UFDs at the redshifts of the Cooke et al. (2014) 
DLAs, and that some were forming their first burst of 
stars, while others were forming their second burst. 

We now compare the DLA population to the stars in 
Leo IV and Canes Venatici H, the two syste ms which are 
likely to form stars at the redshifts of the [Cooke et al.j 
(2014) DLAs. The results are shown in the bottom panel 
ot h'ig. [T^ The 4 stars in Leo IV are consistent with 


the [a/F^ abundances of the DLAs, while of the 6 stars 
in Canes Venatici H with [Fe/H]< —1.5, 4 show [a/Fe] 
enhanced significa ntly above that of th e DLAs. It should 
be noted that while Weisz et al. (2014b) find that the best 
fit for the star forrning hist ory ot Caries Venat ici H is a 


model with multiple bursts. Brown et al. (2014) argue for 
a single burst. We present tne following possible scenario 
for its evolution. The stars with [a/Fe] ^ 0.8 formed 
in the first burst, followed by reionization or supernova 
feedback turning star formation off. The galaxy accreted 
low-metallicity gas which was then polluted by Type la 
supernovae. The second burst of stars therefore formed 
at a lower [a/Fe] for a given [Fe/H]. 

If the DLAs trace gas in UFDs, it is likely that we are 
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- 
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^ 0-5 
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- 4.0 - 3.5 - 3.0 


- 2.5 - 2.0 

[Fe/H] 


- 1.5 - 1.0 


Fig. 11.— [a/Fe] vs [Fe/H] for our simulation f or five simulation ru ns (black points), the [Vargas et al.|(|2013|l sample of stars in UFDs 
in 0.2 dex metallicity bins (red triangles), and the [Cooke et al.]l|2014^ DLAs (blue squares). 


observing some of them during or before their second 
burst of star formation. This provides an explanation as 
to why the DLAs tend to have lower [a/Fe] at low [Fe/H] 
when compared to the UFDs and the average expected 
from Type II supernovae. High [a/Fe] stars formed in the 
first burst of star formation before and during reioniza¬ 
tion, followed by a pause in star formation. By the time 
of the second burst of star formation, sufficient time has 
passed for Type la supernovae to occur, resulting in the 
reduced [a/Fe] seen in the DLAs. This also provides an 
explanation for low [a/Fe], low [Fe/H] stars in systems 
such as Canes Venatici H and Ursa Major 1. 

In the two-burst mo del, the DLAs in the clump at 
[Fe/H] ^ —2.0 in Fig. 12 are forming their first burst 
of stars, while those in tlie clump at [Fe/H] ^ —2.5 are 
forming their second burst. This scenario predicts the 
existence of [a/Fe] = 0.5-0.7 systems at [Fe/H] ~ —2.5, 
which have not yet been observed. 

4.3. Other possibilities 

A third possibility is that the DLAs do not evolve 
from the clump at [Fe/H] ^ —2.5 to the clump at [Fe/H] 
~ —2.0, but that they instead represent syst ems with dif¬ 
ferent masses and star formation histories. ICooke et al.l 
(2014) found that systems with higher velocity widths 
showed higher metallicities. While it is impossible to 
extrapolate directly to halo mass, systems with lower ve¬ 
locity widths will on average have lower halo masses. The 


a/Fe 

a/Fe 


[Fe/H] ^ -2.5 systems may therefore represent systems 
which correspond to UFDs, while the higher metallicity 
systems may be classical dwarf spheroidals. However, 
this does explain why the DLAs are observed with 
^ 0.3, when Type la supernovae could suppress 
below this level. 

Our models in Section 0 showed that systems such 
as M65CH which do not mow star formation can still 
reach the DLA column density limit along some sight¬ 
lines. Some of the DLAs may therefore be systems that 
are not forming stars at the time they are observed. The 
clump at [Fe/H] ~ —2.5 could be systems which formed 
few stars, with most of the enrichment coming from later 
Type la supernovae, while the [Fe/H] ~ — 2.0 systems can 
be explained by the single-burst model. Observations of 
neutron-capture elements such as barium may be able to 
provide support for or rule out this scenario. Neutron- 
capture abundances could also test the possibility that 
the [a/Fe] abundances result fro m a bimodal Type la 
supernova delay-ti me distribution (Mannucci et al.|2006 


Yates et al. 2013), in which approximately 507o ol the 
supernovae occur at t ~ 50 Myr. This could result in an 
early decline in [a/Fe], followed by a period of few Type 
la supernovae, during which [a/Fe] is flat with [Fe/H], 
as is observed in the DLAs. 

An alternative explanation for the variation in [a/Fe] 
is the enrichment of some DLAs by pair-instability super¬ 
novae. Pair-instability supernovae eject a large amount 
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• Model 

I DLA Observations 
iji UFD observations 



- 3.5 


- 2.5 - 2.0 

[Fe/H] 


- 1.5 


- 1.0 


Fig. 12.— Top: Single burst scenario with truncated IMF at the high mass end. Bottom: Two burst scenario with ac cretion of low 
[Fe/H ] gas after the first burst. [o/Fe] vs [Fe/H] for three simulation runs in each plot is represented by black points, the | Vargas et al.| 
( |2013[ l sample of stars in UFDs in 0.2 dex metallicity bins by red triangles, and the [Cooke et ah'] < |2014[ | DLAs by blue squares. 


of iron, resulting in low [a/Fe]. Wise et al. (2012) sug¬ 
gests this as an explanation for the apparent metallicity 
floor of DLAs at [Fe/H] = —3. The only exceptions are 
two DLAs with high [a/Fe] show [Fe/H] < —3, which 
these may have been enriched by an event with more 
usual supernova yields. In future work we will investi¬ 
gate different types of supernovae, such as hypernovae 
and pair-instability supernovae, to determine their effect 
on our higher mass models. 

Finally, we note that factors outside the scope of our 
model could affect our result. For example, in this work 
we assume that all metal species mix into the gas in the 
same way. While different elements may have different 
dust formation time-scales, the gas and dust physics in¬ 
volved is difficult to track. The best work to date on 
this issue has involved considering the dust and molec- 
ular gas phases in post-proce ssing (Krumholz & Gnedin 


2011 Gnedin & Draine||2014|). These complicated issues 
will be considered m future work. However, none of the 
scenarios in this subsection can resolve the discrepancy 
between the UFDs and the DLAs, requiring different en¬ 
richment histories for the two classes of systems 


5. CONGLUSIONS 

Using simulations of star formation and chemical en¬ 
richment in systems with dark matter halo masses of 
10 ^ Mq, we have modelled the chemical abundances of 


very metal-poor DLAs. Our conclusions follow: 


1. For the DLAs in which Gooke et al. 


(2014) were 


able to derive physical quantities, the measured 
sound speeds, temperatures, gas masses and gas 
densities are mostly consistent with gas in our mod¬ 
elled Mvir = 10 ^ Mq halo, although the temper¬ 
ature and sound speed are slightly too low, sug¬ 
gesting the observed DLAs may have slightly larger 
halo masses. 


2. The column densities of the Gooke et al. (2014) 
DLAs can be reproduced by models ot gas m 
a 10 ^ Mq halo which include feedback from the 
preionization and supernova of a massive star. 


3. Multiple supernovae are req uired to explain the 
metallicities of 21 of the 23 Cooke et al. (2014) 
DLAs, except in the case of a pair-instability su¬ 
pernova as in 


Wise et al. (2012). 


4. Our model of star formation and chemical en¬ 
richment can reproduce [a/Fe] for DLAs with 
—2.5 < [Fe/H] < —2.0. While there is significant 
scatter for individual systems, the average UFD 
abundances also agree. 


5. For —3.0 < [Fe/H] < —2.5 there is some ten¬ 
sion between our model and the DLA abundances, 
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with the DLAs showing mean [a/Fe] = 0.25, while 
our model has [a/Fe] ~ 0.5. [a/Fe] is also sup¬ 
pressed for the DLAs compared to the average stel¬ 
lar metallicity in UFDs and compared to the mean 
expected for Type II supernovae for a Kroupa or 
Salpeter IMF. 

6 . One explanation for the abundances of DLAs is a 
scenario with a truncated IMF, or a selection effect 
where only systems with lower mass supernovae re¬ 
tain their gas and are therefore observed as DLAs. 
However, this does not explain the the suppression 
of [a/Fe] compared to abundances of stars in UFDs. 

7. A scenario that can explain the mismatch between 
the DLA and UFD abundances assumes that the 
[Fe/H] = —3.0 to —2.5 DLAs have been enriched 
by two bursts, while the higher [Fe/H] DLAs have 
been enriched by only one. Two-burst DLAs form 
their first burst before reionization before losing 
all their neutral gas. At a later time, they accrete 


low-metallicity gas and commence a second burst of 
star formation. Type la supernovae from first burst 
stars enrich the gas without a delay time. However, 
this scenario does not explain why no DLAs have 
been observed at [Fe/H] ~ —2.5 with high [a/Fe]. 

8 . It is possible that the [Fe/H] = —3.0 to —2.5 and 
the [Fe/H] = —2.5 to —2.0 DLAs have very dif¬ 
ferent star formation histories. For example, the 
low [Fe/H] systems may have formed few stars and 
therefore experienced few or no Type H super¬ 
novae. These systems were then enriched at a later 
time by one or more Type la supernovae, resulting 
in the low observed [a/Fe]. The higher [Fe/H] sys¬ 
tems would then be explained by the single burst 
scenario. 


We thank the anonymous referee for useful comments 
which greatly improved the clarity and quality of this 
work. DW is funded by an Australian Postgraduate 
Award. 
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